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Abstract. I review five of Bolidan Paczynski's ideas on black hole accretion 
disk theory. They formed my understanding of the subject and often guided 
intuition in my research. They are fundamentally profound, rich in physical 
consequences, mathematically elegant and clever, and in addition are useful in 
several technically difficult practical applications. 



1. Introduction 

1.1. Motivation: the five "easy" pieces 

This review article does not give a full (or even coherent) description of black 
hole accretion disk theory. Instead, it only reflects briefly on a few of Bohdan 
Paczynski's particularly important contributions to the subject. Bohdan had 
several brilliant ideas about black hole accretion. Today, most of them are 
mainstream and firmly accepted. One is still considered controversial. 

1. The Potential: For free particles, both Newton's and Einstein's orbital 
dynamics are described by the same principle, and indeed the same equation: the 
orbital frequency follows from the first derivative, and the epicyclic frequencies 
follow from the second derivatives of the "effective potential". Therefore, as 
Bohdan pointed out, by a proper definition of an artificial Newtonian potential, 
one should be able to accurately describe (formally in Newton's theory!) the fully 
general relativistic orbital motion. He then guessed the form of the "Paczynski- 
Wiita" potential which immediately became a great success. The potential is 
simple and very practical in numerous applications. 

2. The Doughnut: One is often interested in phenomena that occur on a 
"dynamical" timescale much shorter than the "viscous" time tc needed for 
angular momentum redistribution, and the "thermal" time ts needed for entropy 
redistribution. In modeling such cases, Bohdan noticed, the angular momentum 
and entropy distributions may be considered free functions. Therefore, when 
studying processes with <C to = Min(t£,t5), one may ad hoc assume angular 
momentum and entropy distributions. A constant angular momentum is the 
simplest possible assumption. It was used by Bohdan and his Warsaw team to 
construct the "Polish Doughnuts", i.e. super-Eddington thick accretion disks. 

3. The Funnel: In Warsaw we found that a long, narrow empty funnel 
forms along the axis of rotation of a Polish doughnut. It collimates radiation to 
hyper-Eddington fluxes. 

4- The Roche Lobe: In equilibrium, the isobaric surfaces of the accretion 
disk coincide with the surfaces of constant effective potential, called the "equipo- 
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tential" surfaces. Near a black hole one of the equipotential surfaces, called the 
"Roche lobe", self-crosses along the "cusp", i.e. a circle r = Vcusp located be- 
tween the innermost stable circular orbit at r = rjsco^ a-nd the marginally bound 
orbit at r = r^b- Bohdan's knowledge of close binaries helped him to realize 
that the black hole Roche lobe overflow must induce dynamical mass loss from 
the disk. Thus, for r < Vcusp accretion is not caused by "stresses" that remove 
angular momentum, but by the strong gravity. The "inner edge" of an accretion 
disk is at the ISCO only for accretion flows with very small M /M^dd that are 
radiatively efficient (Shakura-Sunyaev). He understood that the efficiency of 
accretion rj (defined by L = rjMc^) equals the Keplerian binding energy at the 
location of the cusp: when rcusp = risco, efficiency takes its maximal possible 
value, but it goes to zero when rcusp approaches rmfe- All radiatively inefficient 
accretion flows (RIAFs) like advection-dominated accretion flows (adafs) and 
slim disks have their inner edges close to r^b- 

5. The Inner Torque: Against the view of an influential part of the black 
hole community, Bohdan argued that when the accretion disk is vertically very 
thin, for very small accretion rates the viscous or magnetohydro dynamic (MHD) 
torque at the inner edge should be vanishingly small. 

1.2. Analytic models and computer simulations 

Bohdan was a master of massive supercomputing methods in astrophysics. Some 
of the computer codes that he developed, e.g. for stellar structure and evolution, 
are the standard tools in the field. However, although supercomputer MHD sim- 
ulations are a large part of black hole accretion research today, Bohdan was not 
himself involved in this activity. Indeed, all Bohdan's ideas that I discuss here 
are based on simple analytic models. One may ask — why? This I do not know, 
but I can think of a possible reason. The present-day MHD simulations do not 
address the fundamental problems that the theory of black hole accretion faces. 
Our understanding of black hole accretion rests on analytic models. Indeed, 
all the fundamental features of black hole accretion flows that have been calcu- 
lated in MHD simulations, were previously well understood in terms of analytic 
models. It seems that today, and in the foreseeable future, the super-computer 
simulations will keep confirming rather than solving or discovering. Today, the 
most important message of the supercomputer simulations seems to be that the 
approximations and simplifications adopted in the analytic models have been 
rather wisely chosen (see Figured] to illustrate this point). 

1.3. Notation and units 

In this review, I use the Kerr metric in the standard spherical Boyer-Lindquist 

coordinates t, (f), r, 9 with the H signature. The Kerr metric is stationary 

and axially symmetric. 



i.e. its metric tensor components are (known) functions of radial and polar 
coordinates r, 6, the mass M and the dimensionless spin parameter \a\ < 1, 



ds^ = gtt df + 2gt^ dt d4> + g^,^) d(j? + grr dr"^ + gee dO'^ , 



(1) 




(2) 
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Figure 1. Isobaric surfaces in a black hole accretion flow. Left: according to 
a very simple, analytic Polish doughnut model. Right: according to a state-of- 
art, full 3D MHD numerical simulation. Figure taken from the Living Review 
on black hole accretion disks by Abramowicz & Fragile, in preparation. 
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However, in all specific examples, I will use a much simpler Schwarzschild metric 
that describes spacetime geometry of a non-rotating (a = 0) black hole. 



r 



1-^1 dr 

r 



2 2 

r 



d9^ + sm^ 9 d(b^ 



with the gravitational radius re is defined by, 

2GM 
ro = — 5— 



(3) 



(4) 



My apologies to Bohdan; here I often but not always use the c = 1 = G units 
he greatly disliked. 



2. The Potential 

According to Newton's theory of gravity, the motion of a free particle in a 
stationary, axially symmetric gravitational potential ^{r,9), is confined to a 
plane, which may be chosen as the "equatorial plane", 9 = Tr/2. Circular orbits 
in the equatorial plane are characterized by 

vr 

r = ro = const, 9 = 9q = — = const, (5) 

and the slightly non-circular, slightly off-plane orbits are characterized by 

Sr{t) = r(t) - ro « ro, d9{t) = 9{t) - 9o « 9o. (6) 

The specific (per unit mass) energy £ = [V'^ + (v^)^]/2 -|- <I> and the specific 
angular momentum C = rv^ are constants of motion. Here V"^ = (fr)^ + {vef' = 
(5r)^ -|- {59)"^ <C {v^)^ ■ One defines the effective potential 



2 r^ sin^ 9 

and writes, 



U{T,9,C) = ^r,9) + ,^:j^—^. (7) 



U^'' = £-U{r,9,C). (8) 

The Keplerian angular momentum distribution C = Cxir) in strictly circular 
orbits ([5]) follows from the condition, 



dU 

dr 



= 0, (9) 

c 



while the deviations ([6]) obey the simple harmonic oscillator equatiorQ 

5¥ + ujl5r = 0, 69 + lo159 = 0, (10) 



^When deriving the "radial" equation, one puts V = dr/dt = 5r in and when deriving the 
"vertical" equation one puts V = dO / dt = 59 in (|8]) . 
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with the "radial" ujr, and the "vertical" ujg epicyclic frequencies given by, 

So much for the Newtonian case. The fully relativistic case is remarkably similar. 
The relativistic analog of the Newtonian condition for strictly circular orbits ([5|) 
may be written as a condition for the components of a particle's four-velocity 

= dx^/ds, 

u' = {u\u't'), u'^ = 0, u^ = 0. (12) 
The nearly-circular orbits are characterized by, 

u' = {u\u^, u\ u^), {u'\u'^) < (n*, u't>). (13) 

As in Newton's theory, the energy and angular momentum are two constants 
of the free (i.e. geodesic) motion. Their standard definitions read 8* = ut 
C* = —u^, but for our purpose it is convenient to rescale them according to, 

^Kerr = ln£^*, /^Kerr = (14) 

Obviously, £^Kerr and ^Kerr are also constants of geodesic motion. It is also 
convenient to define the small "deviation" velocity V by its square. 



Vierr = 2 In 



1 - grriun- geeiuy « 1, (15) 



which is obviously always positive because —grr > and —gge > 0. Then, most 
importantly, one introduces the effective potential, 

Z^Kcrr = -\ In (5" - 2£g'^ + £2 gH^ . (le) 

With these definitions, one may easily cast the general equation UiUkg^^ = 1 
into the form 

2^Kerr = ^Kerr — ^Kerr(?', 0, Z^Kerr), (17) 

which is identical with the corresponding Newtonian equation ([8]). All quanti- 
ties that appear in (jl7p have the same physical meaning as the corresponding 
Newtonian quantities, and go to the right limits when v^p/c <C 1, V/c <C 1, 
{2GM/(P')/r = rc/r ^ 1- The same equations have the same solution^ and 
therefore one may just use the same procedure as in Newton's theory to get the 
Keplerian angular momentum from the first derivative of the effective potential, 
i.e. from equation Q, and the epicyclic frequencies from the second derivatives 
of the potential, i.e. from equation (fTT|) . This is indeed the standard, text-book 
way used in general relativity. For example, in the special case of a non-rotating 



^This well-known remark by Richard Feynman is "unsourced" according to wikiquote. 



6 




Figure 2. The Paczynski-Wiita p otential ((2T|) was first introduced in the 
paper bv lPaczvhski fc Wiital ()1980[ ). This is the most frequently cited paper 
by Bohdan on the subject of the black hole accretion. It ranks in the fourth 
place on the list of his most frequently cited papers. The figure and citation 
numbers were taken from the SAO/NASA ADS Astronomy Query service on 
June 6, 2008. 



(i.e. Schwarzschild) black hole, the relativistic effective potential ()16p takes the 
form (on the equatorial plane), 



^Sch 



1 

--In 
2 



1 _ !^ 



-1 



£2 



(18) 



and demanding that its first derivative is zero leads to the expression for the 
Keplerian angular momentum, 



K 



GMr^ 
(r - rcY" 



(19) 



Bohdan noticed that the same expression for the Keplerian angular momentum 
one gets for the effective potential given by the formula (on the equatorial plane) , 



Upw 



GM C?_ 

r — tq 2r2 



(20) 



The second term on the right hand side, C? I2r^^ is the standard Newtonian 
centrifugal part of the effective potential. The first term should be therefore 
identified with a gravitational potential. 



PW 



GM 

r - re' 



(21) 



which is indeed the celebrated "Paczynski-Wiita" or "pseudo -Newtonian" poten- 
tial, f irst introduced by Bohdan in his paper with Paul Wiita (jPaczvnski &: Wiita 
19801 ). 
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3. The Doughnut 



Several interesting phenomena that occur in black hole accretion disks have 
characteristic timescales t^, much shorter than the time tc needed to redistribute 
the angular momentum, or the time needed to redistribute the entropy 



t* <C to = Mm {tc,ts) 
During the time t^:, distributions of angular momentum and entropy, 

£ = £{r,e), S = S{r,e) 



(22) 



(23) 

do not change. The functions (|23|) cannot be unambiguously calculated from 
first principles, as they depend on unknown initial conditions, and on unknown 
details of dissipative processes such as e.g. turbulence, convection, or radiative 
transfer. Therefore, when they are calculated, for example with the help of 
lengthy and difficult numerical simulations, results are arbitrary, at least to 
some degree. 



/ step 1 \ 

arbitrary 
\ assumptions J 



I step 2 \ I step 3 \ 



difficult 
y calculations / 



C = C{r,e) 
\S = S{r,e) J 



(24) 



Paczynski argued that there is no obvious reason to believe that it would be 
easier to guess what needs to be assumed in step 1 than to assume (guess) the 
result in the step 3. In the late 1970s and early 1980s he and his collaborators 
explored accretion disk models in which distributions of angular momentum and 
entropy (p3]) were assumed ad hoc. 

For a perfect-fluid matter, T*^ = (p + ^)u'^ Uk — P^^k^ with p and e being 
the pressure and the total energy density. Equations of mass conservation and 
energy-momentum conservation take the form. 



0, 



(25) 



where p is the rest-mass density. From these equations one deduces that two 
Bernoulli-type quantities are constant along the stream-lines of the perfect fluid 
matter moving in a stationary and axially symmetric spacetime ([I])-©. 



p + e 
P 



ut, L* 



p + e 
P 



(26) 



Thus, the purely kinematic quantity, defined previously as the specific angular 
momentum for particles, 

r = — = 



(27) 



is also constant along fluid's stream-lines. If the fluid moves along circular 
stream-lines ([12]), the equation of motion VjT*^ = yields (c.f. equation (fT6]) ). 



yjP_ ^ l Viff**-2£V,ff*'^ + £2V,g'^'^ ^ f dUKcrA 



(28) 
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which may be transformed into, 



+ ^ (29) 



Here = u'^/n* is the angular velocity, related to the angular momentum £ by 
the following formulae, which are very useful in practical calculations. 

^ _ ^ g(/.(/> + Qt4> ^ _ ^gtt+ gt^ ^^^^ 

3.1. The von Zeipel fluid tori 

Suppose that p = p{e), as it would be for isentropic, polytropic or barytropic 
fluid. It this case, the left hand side of equation ([29]) is a perfect gradient, and 
so must be the right hand side, and this is possible when C = C{il). Similarly, if 
C = C{i}), then the right hand side is a gradient, and so must be the left hand 
side. This statement is known as the von Zeipel condition, 

[p = p{e)]^[C = C{n)], (31) 

or, the p{r, 9) = const surfaces conincide with those of e(r, 6) = const, if and 
only i/the surfaces C{r,9) = const coincide with those Q{r,9) = consfl 

The simplest example of a fluid fulfilling the von Zeipel condition (|3ip is 

C{r,9) = Co = const, S{r,9) = Sq = const. (32) 

In this case, equation (|29p has a trivial first integral 

W{p)=l(KeAr,e), (33) 

where W{p) is the enthalpy, defined for isentropic fluids ([32]) by 

W = W(p) = I (34) 
J p + e 

Note that ([33p gives the location of the isobaric surfaces in terms of an explicit 
analytic function of coordinates r, 9. From (|28p it follows that the pressure 
maximum is located at the circle r = rg, 9 = 9q = 7r/2, given by the condition, 

fdUKerr\ ^ / ^Z^Kcrr \ f^r,. 

Let us express the effective potential near the pressure maximum circle tq, 9o by 

AU = Z^Kcrr(r, 9) - UKeAro, Oo). (36) 



Because the first derivative of the effective potential vanishes at the maximum 
pressure circle, and the second derivatives equal to the epicyclic frequencies (c.f. 
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Figure 3. Contours of equipotential surfaces on the meridional cross section 
of a torus with constant angular momentum, exactly calculated from ([33]), 
according to Newton's theory with the —GM/r gravitational potential. The 
contours approach circles around the locus of the maximum pressure, at the 
r — tq, 6 ^ 6{) ~ t: /2 circle, which agrees with the approximate "slender 
torus" solution given by equation (j37p . 



equation (fTT]) l. one concludes that the shape of a "slender" torus, i.e. with 
Ar = r — ro ^ ro, A6' = 9 — 9o <^ Oq, is given by AU = const, with 



1 r 



Ar 



< 1. 



(37) 



Equation (|37p is valid both in Newton's and Einstein's theory. In Newton's 
theory with —GM/r gravitational potential, one has iVr = cvg. Thus, in this case, 
one concludes from equation (j37p that slender constant-angular-momentum fluid 
tori have circular cross-sections. This add itional symmetry was first noticed by 
Bohdan in his paper with Jurek Madej ( Madei Sz Paczvnski 19771 ). and later 
used by iBlaed ( 19851 ) to calculate normal modes of oscillations of slender tori. 
Blaes demonstrated that there exists a set of non-axisymmetric (i.e. with m ^ 0) 
global modes ~ exp(— zcJmni + m^P) with the eigenfrequencies. 



0"n 



-mVtK 



l + ie 



2n + 2 



n 



1,2,3,, 



(38) 



Because ^{amn 
instability is ~ 



) < these oscillations are unstable. The growth rate of the 
tuVLk^ and thus the instability is a dynamical one. Indeed, this 



■^The best known Newtonian version of the von Zeipel condition states that for a barytropic fluid 
p — both angular velocity and angular momentum are constant on cylinders, il = 
C = C{R), with R = r sin 9 being the distance from the rotation axis. 
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is the famous Papaloizou-Pringle inst ability, discovered in the se minal paper by 
John Papaloizou and James Pringle ( Papaloizou &: Pringld T984I ) . 

In the last twenty years, studies of slender tori oscillations have brought 
a lot of important results. Today, these results are concentrated in non-linear 
excitation, damping, and resonances of global epicyclic modes that are thought 
to be relevant to the 3:2 twin peak QPOs. The velocity pattern for the vertical 
epicyclic m ode of a slender torus in t he Kerr geometry is shown in Figure [U 
taken from lStraub &: Sramkoval ( 20081 ). 




Figure 4. Velocity pattern for the vertical epicyclic mode of oscillation of a 
sle nder torus around a Kerr b lack hole as calculated from an analytic model 
bv lStraub fc Sramkoval (|2008h . 



3.2. General case 



Jaroszyhski &: al ( 1980l ) discuss a general, non-barytropic case of angular mo- 
mentum and entropy distributions (j23p . They start by noting that ii 9 = 9{r) 
is the equation for the isobaric surface p(r, 9) = const, then obviously 



dr 



drp 
dep' 



(39) 



They further note that the right hand side of the above equation may be cal- 
culated from (|28p . if one writes it twice for i = r and i = 9, divides the two 
equations, and changes the sign. 



drP _ dr g** - 2£ drg^"^ + dr-g^ 
~ ~deg*^-2Cdeg^^ + C^ 



(40) 



The point here is that because the Kerr metric g"^^ is known explicitly in terms 
of r and 9, and because the distribution of the angular momentum is given by 
the assumption (|23p . the right hand side of (j40p is an explicitly known function 
of r and 9. One may denote this function by F{r,9), and combine equations 
(j39p and (j40p into the standard form of an ordinary differential equation. 



f = F{r,9), 
dr 



(41) 



with the explicitly known right hand side. It may be directly integrated to get 
all possible locations of of the isobaric surfaces. 
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Figure 5. Shapes of Polish doughnuts. Upper row: for ahnost constant 
angular momentum, /3 w 0. Lower row: for almost Keplerian angular mo- 
mentum, /3 w 1. The parameter 7 increases in equal steps from left, 7 = 0, to 
right, 7=1. 



Qian &: al ( 20081 ) have recently argued that a physically flexible and realistic 



angular momentum distribution may be assumed of the form, 



£(r,( 



£o(^fP)^ (sin0)2T for r > 
Cms (sin 6*)^^ for r < 



(42) 



where Lms = [Lxifms) / i^a]^ ■, and Cq, P, and 7 are the three free constant 
parameters of the distribution. Figure [5] shows several models calculated for 
different choi ces of £0, /3, 7. 

Recentlv lKomissarov ( 20061 ) constructed models of magnetized Polish dough- 
nuts. 



4. The Funnel 

Long, narrow funnels along the axis are a genuine feature of Polish doughnuts. 
Bohdan realized that they may be relevant for beaming radiation to highly 
super-Eddington fluxes and collimating relativistic jets. 

4.1. Super-Eddington luminosities in funnels 

In this section, I will use Newton's theory and cylindrical coordinates r, z, (j). In 
Newton's theory, the first integral (j33p that describes equilibrium of the constant 
angular momentum, constant entropy fluid takes the form, 

GM 1£2 

- + -— + W{p) = const. (43) 

(r2 + z2)2 2r 

This equation cannot be obeyed at the rotation axis (where z 7^ 0, r = 0), which 
has an obvious, but important, consequence: no isobaric surface may cross the 
axis. For a constant angular momentum fluid, isobaric surfaces must be toroidal, 
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or open. The marginally open surface has just one point (r = oo,z = 0) at 
infinity. This particular surface encloses the largest possible torus. From (j43p it 
is obvious that in this case W{p) = W{0) = 0. The maximal pressure is located 
at a circle z = 0,r = vq = GM/Cq. Using the radius vq as a scale, ^ = r/ro, 
r] = z/vq, w = W/{GM/rQ), one may write equation (|43p in the dimensionless 
form, and solve for = r/(^, w) to obtain, 

^■^ = Q\^) = Ai^(i-2ewy^ -e, -i/2<w<ws. m 

The value w = —1/2 gives the location of the center, and w = ws < the 
location of the surface. For the slender torus ws ~ —1/2, and the fat torus 
Ws ~ 0. Very far from the axis of rotation one has \2w(,\ ^ 1. Inserting this 



Figure 6. Contours of equipotential surfaces on the meridional cross section 
of a fat torus with constant angular momentum, i.e. with Wg ~ 0. Close to 
axis the isobaric surfaces form a pair of conical funnels. 

into ()44p one gets 7?^+^^ = which means that far from the rotation axis the 
equipotential surfaces are spheres with radius l/jit'l, and that the outer radius 
of the torus is at rout = ro/\ws\- Spherical equipotentials are in accord with the 
fact that very far from the axis, ^ ^ 1, the centrifugal force ~ Cq/^^ is negligible 
with respect to the gravitational force ~ GMj^ . Therefore, the effective gravity 
is determined by Newton's attraction alone, as for spherical stars. For the same 
reason, the radiation power from this part of the surface of a fat, radiation- 
pressure-supported torus (i.e. Polish doughnut) is one Eddington luminositv, 
the same as from a spherical non-rotating, radiation-pressure-supported stao 



*Martin Rees told me as early as 1980 that to him this implied that Polish doughnuts, with 
at least one Eddington luminosity from whatever direction, were obviously too luminous to 
describe very "dim" active galactic nuclei, such as radio galaxies: "while apparently supplying 
tremendous power to their extended radio-emitting regions, the nuclei of most radio galaxies 
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Note, however, that the asymptotically spherical shape of a fat torus is a direct 
consequence of the assumption £(r, z) = const, which was made ad hoc. If one 
adopts a more physically realistic assumption that asymptotically C{r, z) = Ck, 
as in our model (j42p . one may use the standard Shakura-Sunyaev model in its 
radiation pressure version to get the asymptotic shape of the fat torus, 

M. (45) 



Sircmp 



Closer to the axis, \2w(,\ ~ 1, which means that ~ —l/2w, and this to- 
gether with ([li]) gives t/^/^^ = (l/|'w|) — 1 = (1/sin^^) — 1, i.e. that closer 
to the axis, equipotential surfaces corresponding to u) ~ ~ have con- 
ical shapes with the half opening angle ~ ^/w. The surfaces are highly 
non-spherical because centrifugal force dominates. Int egrating effective gray - 
ity along the conical funnel is elementary, and one gets ( Jaroszvhski fc: alll98Cll ) 



that Lg'j^/LEdd = (1/2) ln(l/|ii;5|). This estimate may be used to find the total 
luminosity for the radiation-pressure-supported fat torus, i.e. a Polish dough- 
nut, 

=l.l5log(^-^]. (46) 
-^^Edd ^Edd 2 V ^'o / \ ro J 

It should be clear from our derivation that the logarithmic scaling of the lumi- 
nosity with the torus size is a genuine property of Polish doughnuts, including 
those that have an angular momentum distribution that is not constant. The 
logarithm in (|46p is of crucial importance, as it prevents astrophysically realistic 
doughnuts (i.e. with Tout/'^o < 10^, say) from having highly super-Eddington 
luminosities. Thus, the theory predicts that for such "realistic" fat tori, only 
slightly super-Eddington total (isotropic) luminosities, Ltotai ^ 7LEdd; may be 
expected. 

However, because the funnels have solid angles ~ ro/rout) radiation in 
the funnels may be, in principle, collimated to highly super-Eddington values 
-^couZ-^^Edd = ~ ro^t/vQ » 1. This simple estimate agrees w ith a in ore detailed 
modelling of the radiation field of the Polish doughnuts by ISikora ([198 li ) and 



Madaul (|l988l ) who obtained 6 < 10^ for disks with rout/^^o ~ 10^. A typical 



value that follows from observational estimate s for non-blazar active galactic 
nuclei (e.g. iMalkanI Il989l : ICzernv &: Elisl \l98l\ ) is 6 ~ 10, but of course for 



blazars and other similar sourc es, e.g. for ULXs, if they are powered by stellar 
mass black holes, as arg ued bv lKingl (|2002l 1. it must be 6 > 10. 



Such high v alues of are also coi isistent with the id ea, suggested by 
Paczvnski ( 1980l ) and independently by iLvnden-Bell ( 19821 ). that relativistic 



electron-positron e e"*" jets may be very effectively accelerated by the radiation 



pressure in the fat tori funnelqj. Note that if the flux in the funnel is times 



emit little detectable radiation." He and his collaborators at Camb ridRe later fou nd a possible 
solution to this puzzle in terms of the ion-pressure-supported tori Re es fc all (|l982r ). The ion tori 
have the same shapes as Polish doughnuts (in particular funnels) but have much lower, indeed 
very sub-Eddington, luminosities. The power in jets comes from tap ping the rotati o nal en ergy 
of the central black hole by the Blandford-Znajek mechanism iBlandford fc Znaiekl l|1977f ). 



^Lynden-Bell called this an "entropy fountain". 
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the Eddington flux, the e~e+ plasma feels the "effective" radiative force corre- 
sponding to the Eddington ratio mp/mp, ^ 10^ times greater . Detailed calcula- 
tions (e.g. Abramowicz &: Piranlll980l : Sikora Sz Wilson 1981 : Abramowicz &: al 



Il990l ) demonstrate that indeed the e~e+ jets may be accelerated in funnels 
up to the Lorentz factor 7 < 5. However, if jets are initially pre-accelerated 
by some bl ack-hole electrodynamical processes (such as the Blandford-Znajek 
mechanism, Blandford Sz Znaiek 19771 ) to highly relativistic velocities 7 > 10^, 



they will be decelera ted in the funnels by C ompton drag, reaching the asymp- 
totic Lorentz factor ( Abramowicz &: aj||l99Cll ). 



1 

m„ \ 3 _ 1 



=ioxei (47) 



Observations show that 700 < 10 , and thus equation ()47|) suggests that Q < 10 
4.2. The physical reason for super-Eddington luminosities 



I will n ow review a very general argument, first presented bv l Abramowicz. Calvani Sz Nobili 



( 1980l ). which shows that the upper limit for the luminosity of rotating bodies 



in equilibrium must differ from the Eddington limit. I will consider two "astro- 
physical" cases: 

Rotating stars. The surface of the star has the topology of a sphere. The 
whole mass M is included in the sphere. 

Accretion disks. The surface of the disk has the topology of a torus. The 
mass of the disk Mdisk is contained in the torus, but the mass Mcentr of the 
central black hole is outside. The total mass M = Mcentr + M^isk ~ -Mcentr, 
because Mcentr ^ Mciisk- In accretion disk theory it is customary to neglect the 
mass of the disk, so formally M = Mcentr, and Mjisk = 0. The Eddington limit 
always refers to the total mass M. 

Let 

frad t)6 the local flux of radiation somewhere at the surface of the 
body, S. The corresponding radiative force is Fj-ad = ('7rad/c)frad- Let Feg = 
Fgrav + Froti be the effective gravitational force, with Fgrav = mV^ being the 
gravitational force ($ is the gravitational potential), and with Fj-ot = m{Q'^r) e^ 
being the centrifugal force (O is the angular velocity, r is the distance from the 
axis of rotation, and e,- = Vr a unit vector in the off-axis direction). 

The necessary condition for equilibrium is Fiad < Fcff. Prom this one 
deduces the Eddington limit for rotating perfect-fluid bodies, 

L = / f,ad • dS = — / Frad • dS < — / Feff • dS = Lg-^^d (48) 
Js 0-rad Js 0-rad Js 

Using Gauss's theorem to transform the surface integral of Feff into a volume 
integral of V-Feff , Poisson's equation V^^> = AnGp , and introducing the speciflc 
angular momentum £ = J7r^, one gets after twenty or so lines of simple algebra. 



-^Edd — ^Edd 



-^mass ' shear -^vorticity 



(49) 



where LEdd is the standard Eddington limit for a non-rotating star. 



LEdd = 1-4 X 1038 ( ^ ) [ erg sec"!] . (50) 
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and where the dimensionless number -^^ass depends on whether the body is a 
star, or an accretion disk, 



-, » ( 1 for stars, 

Xl^s = T7 i_PdV={ }. (51) 

for accretion disks. 



M 



V 



Here ^s\ear -'^vorticity dimensionless, necessarily positive quantities, equal 
to shear and vorticity integrated over the whole volume of the body, 

x;„,.., = :f^/^r-(V£.V£)dF. (53) 

Shear increases the Eddington limit, and vorticity decreases it. 

The rotation of astrophysical objects is far from simple, but insight can 
be gained by considering a simple power law for the angular momentum dis- 
tribution, £.{r, z) = £o^"; with Cq and o constant. Rigid rotation has a = 2, 
Keplerian rotation a = 1/2, and constant angular momentum rotation a = 0. 

vorticity 



-^shear Emd --^vorticity 

are related to a by -'^shear /^vorticity = - ^Y/a?- This 



means that X^^^^^ > X^^rticity when a < 3/2. 

Rotating stars have ^g^ear < ^vorticity > because they rotate almost rigidly. 
Thus, for rotating, radiation pressure supported stars, L ^Edd < -^Edd always. 
Contrary to this, tori with constant angular momentum are dominated by shear, 
"'^s^hear ^ 1 ^ ^vorticity ^^^^ Consequently, when they are radiation pressure 
supported, L w U^^^ > L^dd- 

4.3. Rise and Fall of the Polish Doughnuts 

At the time of their debut, Polish doughnuts could theoretically confirm the ob- 
served super-Eddington luminosities, highly collimated beams of radiation, and 
perhaps even the relativistic speeds of jets, fulfilling the principle attributed 
to Eddington: one should never believe any experiment until it has been con- 
firmed by theory. These virtues initially attracted some interest in the astro- 
physical community, but the interest quickly drained with the discovery of the 
Papaloizou-Pringle instability. It was thought that Polish doughnuts must nec- 
essarily suffer from the instability and thus, in reality, they cannot exist. The 
important discovery bv lBlaej ( 19871 ) that the Roche lobe overflow stabilizes Pol- 



ish doughnuts against the Papaloizou-Pringle instability came too late — in the 
advent of numerical simulations of black hole accretion flows. Too late, because 
numerical simulations rediscovered and absorbed many of the results of Polish 
doughnuts. Today, these results exist in the consciousness of many astrophysi- 
cists as a set of several numerically established, important but unrelated facts. 
They do not form a consistent scheme that the Polish doughnuts once offered: 
clear, simple, following directly from the black hole physic^ 



J6^ 



® There are at least three brilliant an d very pedagogic a l exp ositions of the Polish doughnuts 
scheme: two by Paczyhsk i himself (|Paczvhskil 1 19821 . 119981 ) and one in the text book by 
iFrank. King fc Raind (20021 ). 
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5. The Roche Lobe 



A powerful analytic model for accretion flows near the ISCO was worked out 
in terms of the flow equipotential st ructure by Paczynski a n d collaborators 
thirty years ago in Warsaw (see e.g. Abramowicz Sz all Il978l : iKozlowski Sz all 



19781 : 1.Taroszvhski fc allll98nl : lAb ramowic zjll98lll. The Warsaw model accurately 



describes the flow hydrodynamics for the standard Shakura-Sunyaev thin disks, 
slim disks, adafs, ion tori, and thick Polish doughnuts. It impressively agrees 
with m ore recent numerical simulations of accretion flows, as directly checked 
e.g. by llgumenshchev Sz Abramowicd ( 20001 ). 

In an equilibrium described by the Bernoulli equation ()43p . surfaces of con- 
stant enthalpy, pressure and density coincide with surfaces of constant effective 
potential, U{r,z) = const. The surface of the disk is given by P = const = 0. 
Thus, equilibria may only exist if the disk surface corresponds to one of the 
equipotentials inside the Roche lobe, i.e. in the region indicated by yellow in 
Figure |71 If the fluid distribution overflows the Roche lobe, i.e. the surface of 
the disk for r » ri coincides with U{r,z) = Us < Z^/q,, the equilibrium in the 
region r < ri is not possible, and the disk will suffer a dynamical mass loss, with 
accretion rate M\n. 
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5.1. The stationary Roche lobe overflow 

An an alytic formula for M for stationary flows was first calculated bv lKozlowski &: aJ 
( 19781). who used E instein's theory. Here we review another derivation, done by 
Abramowicj ( 19851 ). as we later use the same assumptions and notation to cal- 



culate M for a non-stationary (oscillating) disk. In particular, we assume that 
the Roche lobe overflow is small (quadratic in disk thickness H), 



h{ri,z) 



1 



2 2 
K Z , 



H 



(54) 



that the equation of state is polytropic, and that the radial velocity v"^ connected 
to the mass loss through the cusp equals the sound speed Cg, 



(55) 



In the equations above h* = h{ri, 0) denotes a maximal value of the enthalpy on 
the cylinder r = ri and the subscript L stands for the evaluation of the derivative 
at the point [ri,0]. The local mass flux rii = pv"^ = /i"+^/^/K"(l +n)"n^/^, 
vertically integrated through the cusp thickness and azimuthally around, gives 
the desired total mass flux in terms of the enthalpy, 



M 



2lT 



ridip 



+H 



rhdz 



-H 



(2vr) 



3/2 n 



Kn 



1/2 



K{n + 1) 



r(n + 3/2) +1 
Tin + 2) ■ 



(56) 



Here r(x) is the Euler gamma function. From + h+lA = Us one gets 
AU = Us-U = {l + l/2n) h and from 



dz"^ 



n 



n + 1/2 



dz"^ 



we recover the result obtained by I Abramowicz ( 19851 ). 

M = A(n)— AZ^"+\ 

n+1/2 r(n + 3/2) 



A{n) = (27r)3/2 



1 



K{n+1) 



1 



n+1/2 



r(n + 2) 



(57) 



(58) 



(59) 



Figure [5.1.1 shows an excellent agreement between the analytic formula (j58p and 
results of large-scale, 3D, non-stationary numerical simulations. 

5.2. Non-stationary Roche lobe overflow 

We will calculate M for non-stationary oscillating flows in a special but im- 
portant case, assuming that the poloidal part of velocity v = {v'^,v^) may be 
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Figure 8. The line is according to the analytic formula ([55)1 and 
the points are from non-statio nary, 3D numerical simulations by 
llgumenshchev fc Abramowicj ()200Cl[ ). 
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derived from a potential. In this case, the non-stationary version of the Bernoulh 
equation has the form, 



-^ + — + h + UUs, with v = Vx- 
ot 2 



(60) 



We assume that the osciUation is a smah, non-stationary perturbation to the 
stationary flow considered in the previous section. 



X{r,t) = X{0){r) + ex(i)(.f,t), 



(61) 



where the subscript (0) refers to the stationary flow and the dimensionless pa- 
rameter e <C 1 characterizes the strength of the perturbation. From the defini- 
tion V = Vx one derives. 



vfo) + 2e^(0) • ^^(1) + e'^vfi) 



'9X(i)\^, fdx(i) 



dr 



+ 



V dz 



The perturbed enthalpy profile can be approximated by an expansion 



h = h 



(0) 



(62) 



(63) 



By substituting equations ([62]) and ([63]) into the Bernoulli equation ([60]) and 
equating coefficients of the same powers of e, we get 



1(1) 



1(2) 



h \ dx(i} 




(64) 
(65) 



This way, we express perturbations of fluid quantities in terms of a perturbation 
X(i). To progress further, one must know the function X(i) = X(i)(*;^i-2) that 
describes oscillations. Finding it in general is a difficult global problem. We do 
not attempt to solve it here. Instead, we describe oscillations by an ansatz. 



X{1) = ZVz COSUJt, 



const. 



(66) 



which models a vertical "epicyclic" oscillation with frequency u that rigidly 
moves the fluid up and down across the equatorial symmetry plane. Such a 
rigid mode, with the eigenfrequency equa l to the epicyclic v ertica l frequency, 
io = iOz, has been recently identified by lAbramowicz all (|2006l 1 as a com- 
bination of eigenmodes of slender torus oscillations found previously by iBlaed 
(|1985I 1. This classic japer gives all possible modes, i.e. their eigenfrequencies 
and eigenfuncions, in terms of exact, simple, and explicit analytic formulae. 

The quantity evz = const can be interpreted as the amplitude of the vertical 
velocity. Equations ([M]) and ([65]) give 



= zVzU) sinojt, /i(2) 



--v1 cos^ U)t. 
2 ^ 



(67) 
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The vertical profile of the enthalpy at r = ri is 



h{ri,z,t) = h* 



2 2 

K z + ezVzU) smujt 



1 



-e^vf cos^ iot + Oie'^ 



which is quadratic in the variable z (see the left panel of Fig. [9|) . The position 
of the enthalpy maximum on the cylinder r = ri is shifted from 2; = to height 
6z{t), given as 



6z{t) = 6Zsmujt, 6Z = e 



UJV. 



We can interpret 6Z as the amplitude of the oscillations. Also the value of 
enthalpy in the maximum differs from the stationary case by 



6h* 



h{ri,6z) - h* 



1 2 

-K 

2 



5z' 



(70) 



According to equation (j56|) the instantaneous accretion rate depends on the 
maximal enthalpy as M oc (/i*)"+^. This relation can also be applied in the case 
of vertical oscillations because the z-dependence of enthalpy on the cylinder 
r = r\ remains quadratic and the oscillations do not contribute to the radial 
velocity of accreted matter. Hence, using equations ([5^ . ([57|) and ([70|) and 
assuming that the frequency of oscillations equals the local vertical epicyclic 
frequency, a; = w^, we arrive at our final result 



bM 



M, 



(0) 



(n + 1 



P 



6h* 



2p 



(1+p) 



6z^_ 6Z^ 



(71) 



where 5M = M — M(o) and p = n/{n + 1/2). The quadratic dependence of 

6M on perturbation implies that the frequency of modulation of M is twice the 
frequency of the disk oscillation (see the right panel of Fig. [9]) . 



6. The Inner Edge 



This is a ve ry controversial point of & fundamental importance. In Bohdan's 
own words ( Afshordi &: Paczviiski 20031 ). "Theory of accretion disks is several 
decades old. With time ever more sophisticated and more diverse models of 
accretion onto black holes have been introduced. However, when it comes to 
modeling disk spectra, conventional steady state, geometrically thin-disk models 
are still used, adopting the classical ''no torq ue " inner bound ary condition at 
the marginally stable orbit at the ISCO.{e.g. iBlaes fclil I2OOII ). Recently, the 
no torqu e condition for geometrically thin disks has been challenged by several 
authors (|Kroliklll999l : lGammiel[T999l : lAgol k Krolikll2000l 1. I did not agree with 
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Figure 9. Left: The vertical profiles of the enthalpy h{ri, z) on the cylinder 
r = ri during vertical disk oscillations (thin lines). The amplitude of oscil- 
lations is SZ = 0.2H and we chose the polytropic index of the fluid n = 1.5. 
The figure captures profiles with enthalpy maxima at Sz ~ 0, 0.IH and 0.2H. 
The enthalpy profile for the unperturbed stationary disk is also shown (thick 
line). Right: The modulated accretion rate from the oscillating disk (thin 
solid line). The accretion rate for the stationary disk is plotted by the thick 
line. Time is rescaled by the period of oscillations. For reference we plot also 
the phase of disk oscillations (dotted line) The accretion rate is modulated 
with twice the frequency of oscillations. 



their claim and presented simple arguments why the no torqu e inner bou ndary 
condition is natural if an accretion disk is geometrically thin ( Paczyiiskil [20 00) . 



but the re feree could not be convinced. Thanks to the electronic preprint server, 
the paper Paczviiskil ( 2000l ) is readily accessible to all interested readers, who 



can judge its validity." 

A more detailed qua ntitative analysis of the inner boundary condition, re- 
fining arguments given by 'Paczynski' ('20001) and making them more pr ecise, was 
done byKfs hordi fc Pa czyilski (2003). They found that Krolik ( 19991') . iGammid 
(|1999I 1. and lAgol k KrolikI (|2000l ) were qualitatively correct: a classical thin 



steady-state disk must have some torque at the ISCO, but the effect is not as 
strong as claimed by Krolik, Gammie and Agol. 
Bohdan's line of argument is as follows: 

• Today, the only models that make solid quantitative predictions about the 
disk spectra that could be compared with observations are the conventional 
steady state, geometrically thin-disk models based on the classical "no- 
torque" inner boundary condition at the marginally stable orbit at the 
ISCO. 

• It is not known today whether steady-state thin-disk type accretion is pos- 
sible at all. For serious technical reasons, not even the most sophisticated 
3D MHD numerical models of accretion can describe very thin flows, or 
include radiative cooling. 

• Therefore, all arguments today concerning the behavior of the steady-state 
flows near the ISCO are incomplete. They all make simplifying assump- 
tions. 



22 



• The most fundamental property of the steady-state flows near the ISCO is 
that they must necessarily be transonic (with the speed of sound properly 
defined to include magnetic effects). The sonic point is the critical point 
in the mathematical sense. It is known from the theory of differential 
equations that certain regularity conditions must by obeyed at a critical 
point. 

• Krolik, Agol and Gammie do not properly consider the regularity condi- 
tions at the critical point(s). This is exactly the reason why the stress they 
calculate at the ISCO is so large. In all models which take a proper care 
of this important piece of mathematics and calculate solutions which pass 
the critical points smoothly, the stress at the ISCO is non-zero, but small. 

A very important num erical work t h at su pports Paczyhski's point of view 
was completed recently by IShafee Sz all ( 20081 ) . It describes results of three- 



dimensional general relativistic magnetohydro dynamical simulations of geomet- 
rically thin accretion disk around a non-spinning black hole. The disk has a 
relative thickness h/r ~ 0.005 — 0.1 over the radius range (2 — 20)GM/c^. In 
steady state, the specific angular momentum profile of the inflowing magnetized 
gas deviates by less than 2% from that of the standard thin disk model with the 
zero-inner-torque assumed. In addition, the magnetic torque at ISCO is only 
~2% of the inward flux of the angular momentum at that radius. Both results 
indicate that magnetic coupling across the inner edge in relatively unimportant 
for geometrically thin disks around non-spinning black holes, which is in accor- 
dance with Paczyhski's ideas. However, until a mathematically correct analytic 
model describing thin MHD accretion fiows near the ISCO becomes available, 
the controversy is likely to continue. 
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